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ABSTRACT 
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1.  Introduction 


Approximate  theories  for  calculating  sound  radiation  from  and  sound  scattering  by  a  simply 
supported  cylindrical  shell,  which  is  contained  between  cylindrical  rigid  baffles,  are  now 
basic  to  acoustics,  see,  for  example,  the  standard  text  [1]  on  fluid-structure  interaction  and 
the  reports  by  James  [2-3].  A  report  by  Jones  [4]  gives  the  theory  of  a  simply  supported 
shell  which  is  covered  by  a  coating  modelled  as  flui'’  Many  assumptions  have  been  made 
in  order  to  render  tractable  the  mathematics  of  the  simply  support^  shell  vibrating  in  an 
acoustic  fluid,  among  which  is  the  assertion  that  modal  coupling  of  the  in- vacuo  structural 
modes  by  the  surrounding  fluid  can  be  neglected,  and  for  certain  modes,  structural  loss 
factors  are  greater  than  radiation  loss  factors.  Thus,  while  recognizing  that  no  more  than 
qualitative  numerical  values  or  physics  can  be  anticipated  from  these  models,  their  primary 
role  is  to  help  contribute  valuable  insight  into  the  physical  mechanisms  that  control  the 
acoustics  of  finite  shells. 

An  increasing  use  of  composite  materials  makes  desirable  analytical  methods  for  predicting 
their  acoustic  characteristics.  Recent  work  by  Skelton  &  James  [5]  has  provided  a  substantial 
theoretical  basis  for  calculating  the  acoustics  of  infinite  cylindrically  layered  media 
comprising  anisotropic  solids  and  acoustic  fluids.  A  semi-analytical  approach  is  used;  axial 
and  circumferential  variations  are  represented  by  Fourier  transforms  and  series,  respectively; 
through  thickness  variations  by  the  finite  element  method.  The  aforementioned  work  is  used 
as  the  foundation  of  the  theory  which  is  given  here  with  a  minimum  of  explanation.  The 
geometry  is  shown  in  Figure  1;  a  simply  supported  cylinder,  comprising  layers  of  anisotropic 
solid  and  acoustic  fluid,  is  contained  between  rigid  baffles.  The  objective  is  to  give  an 
approximate  theory  which  can  be  used  to  predict  the  far  field  sound  radiation  when  the 
excitation  is  a  mechanical  force  and  the  target  strength  when  the  excitation  is  a  plane  wave 
at  oblique  incidence. 

The  work  described  herein  was  done  under  item  AS02EH14  of  the  Strategic  Research 
Programme,  "acoustics  of  novel  materials  and  hull  configurations”,  it  contributes  to  the  third 
Milestone. 
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2. 


Mathematics 


2.1  Fourier  Expansions 

Letting  denote  axial,  tangential  and  radial  displacements,  then 

the  following  expansions  satisfy  the  simply  supported  boundary  conditions  at  the  ends  of  the 
cylinder,  z=-L  and  z=L, 

=  52  «j(r,/i,«)cos[mii(z+Z,)/2Z,]exp(m4)), 

=  52  «4('‘.n.w)sin[mit(z+L)/2I]exp(in(J>),  (2.1.1) 

=  51  tt,(r,/i,w)sinIinit(z+L)/2Z.]exp(i/i<|)), 

where  the  sums  are  understood  to  extend  over  values  of  m  ranging  from  1  to  «,  and  values 
of  n  ranging  from  -«>  to  +«.  The  time  variation  factor  expf-iQr)  is  suppressed  in  these 
and  subsequent  equations.  The  spectral  quantities  are  given  by  the  inverse  transforms,  for 
example, 

♦t  2* 

ujirji,m)  =  (l/2nL)  j  j  t<^(r,(t>,z)sin[mTi(z+I)/2L]exp(-tn4>)d<})dz.  (2  1.2) 

-L  0 


2.2  Dynamic  Stiffness  of  Elastic  Element 

The  procedure  closely  follows  previous  work  [5],  except  that  a  Fourier  series  represents  axial 
variations  rather  than  a  Fourier  integral.  The  element  has  three  nodes  through  its  thickness, 
the  first  coincides  with  the  outer  surface  r=b,  the  second  node  is  at  the  mid-surface 
r={a+b)l2 ,  and  the  third  node  is  at  the  inner  surface  r-a .  Omitting  details,  the  9x9  spectral 
dynamic  stiffness  matrix,  'vhen  the  variations  of  the  displacements  in  the  radial  coordinate 
are  approximated  by  quadratics,  is  found  to  be 

[S(n,m)-<d^M]u(n,m)  =  F(n,m),  (2-2.1) 

where  h 

S(n/R)  =  (A'‘)*^j'B*(r,/i,/n)^DB(r,n,m)rdrA'\ 

(2.2.2) 

* 

M  =  (A-‘)*Y  pN*(r)^(r)nirA-‘, 
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with  the  asterisk  denoting  the  operation  of  complex  conjugate.  The  integrals  are  evaluated 
numerically  by  using  a  two-point  Gaussian  quadrature  formula. 


The  9x1  displacement  and  excitation  matrices,  which  are  the  double  Fourier  senes 
coefficients  of  the  nodal  displacements  and  nodal  stress  excitations,  are 

(2.2.3) 

ujA,njn),u^{a,n/n)M,ia,n,m)Y, 

and 

F(n,m)  = 

Fj(c,n,m),F^(c,n,m),F,(c,n,m),  (2.2.4) 

F.(a,rt,m)^^(a,n,m),F/a,/i,m)]^ 

with  c=(a+b)/2.  The  stress  excitations  are  defined  as  positive  when  they  act  in  the  positive 
directions  of  the  coordinate  axes.  The  9x9  matrix  connecting  spectral  displacements  and 
unknowns  is 

[l  6  000000 

000  I  b  b^  0  0  0 

000000  lb  b^ 

Icc^  000000 

A=00  0  Icc^OO  0 

OOOOCO  Icc^ 

la  000000  (2.2.5) 

000  la  000 

00  0  00  0  1  a  , 
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The  6x9  matrix  connecting  spectral  strains  and  unLiOwns  is 

-o  -or  -or^  0  0  0  0  0  0 

0  0  0  w/r  in  inr  I/*-  1  r 

0  0  0  0  0  0  0  1  2r 

infr  in  inr  a  or  or*  0  0  0 

0  0  0  -1/r  0  r  injr  in  inr 

0  1  2r  0  0  0  o  or  or* 

were  a^mnfZL.  This  matrix  differs  from  the  corresponding  matrix  of  the  infinite  cylinder 
[5],  as  it  reflects  a  Fourier  series  expansion  in  the  axial  direction  rainer  than  a  Fourier 
integral  transform.  The  3x9  polyomial  matrix  is 

1  r  r^  0  0  0  0  0  0 

N(r)  =  00  0  1  r  r*  00  0  (2.2.7) 

0  0  0  0  0  0  1  r  r^ 

The  inverse  of  the  6x6  constitutive  matrix  when  the  material  is  orthotropic  is 

-^rJK  0  0  0 

0  0  0 

-V  /£  -V  l/£  0  0  0  ^ 

j)-I^  ri'rr  n|/rr  '  rr  (2.2.8) 

0  0  0  1/G^^  0  0 

0  0  0  0  1/G^  0 

0  0  0  0  0  1/gJ 

where  the  constants  have  their  customary  meanings  and 

(2.2.9) 

For  fibre  reinforcement  at  a  helical  winding  angle  alor,i  the  cylindrical  surface,  the  matrix 
transformations  to  convert  a  consitutive  matrix  in  a  local  axes  set,  in  which  it  is  orthotropic, 
to  a  global  axes  can  be  found  elsewhere  [6].  Damping  is  simulated  by  allowing  the  elastic 
constants  to  be  complex,  viz  £2£(l-iTi)  etc  where  rj  is  the  hysteretic  loss  factor. 
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After  taking  the  inverse  of  the  9x9  dynamic  stiffness  matrix  in  equation  (2.2.1),  deleting 
rows  and  columns  4,5,6  and  then  inverting  the  resulting  receptance  matrix,  the  6x6  dynamic 
stiffness  matrix  relating  spectral  displacements  and  excitations  at  the  outer  and  inner 
boundanes,  r=b  and  r=a ,  is  obtained  for  the  element  as 

[S,(n^)\ujin/n)  =  F,(n,m),  (2.2.10) 

where 

(2.2.11) 

=  [F^(bjt,m)J^^(b,n,m)J^^(b,n,m), 

are  the  surface  spectral  displacements  and  excitations. 


2,3  Dynamic  Stiffness  of  Fluid  Element 

For  an  acoustic  fluid  layer  of  density  p ,  sound  speed  c ,  inner  radius  a  and  outer  radius 
b,  the  2x2  dynamic  stiffness  matrix  is  obtained,  after  expanding  the  pressure  in  a  Fourier 
sine  series  as 


lilnyb 

ujib^^) 

FJib/t,m) 

yW 

li/nya 

'^|.|(Ya)ff',,|(Yi>) 

«/a,rt,m) 

F^{a/i,m) 

where 


W  =  J',„|(Y6)f/',,,(Ya)  -  J'|„|(Yfl)f/ |„|(y6),  (2.3.2) 

and  ‘f  =  + with  k=<t)lc. 
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For  an  exterior  acoustic  fluid  of  density  p ,  sound  speed  c ,  and  outer  cylinder  radius  fc ,  the  1  x  1 
dynamic  stiffness  matrix  is  obtained  as 


(2.3.3) 


For  an  interior  acoustic  fluid  of  density  p ,  sound  speed  c ,  and  inner  cylinder  radius  a ,  the  1  xl 
dynamic  stiffness  matrix  is  obtained  as 

(2.3.4) 

-p<o^[J|^|(yn)/Y7'l^|(Yn)]u,(n,n^)  =  F^(aji,m). 


2.4  Assembly  of  Elements 

Assembly  of  M  element  matrices,  together  with  matrices  for  the  exterior  and  interior  fluids, 
is  a  straightforward  finite  element  procedure  which  results  in  the  system  matrix  equation 

Z(rt,m)U(n,in)  =  E(n^)  (2.4.1) 

where  Z  is  a  3(Af+l)x3(Af+l)  matrix  of  bandwidth  11.  U  is  a  column  vector  of  the 
3(M+1)  spectral  displacements  at  the  interfaces,  and  E  is  a  column  vector  of  the  spectral 
excitations  which  are  positive  when  acting  in  the  positive  directions  of  the  coordinate  axes. 
A  2x2  fluid  dynamic  stiffness  matrix  can  be  assembled  by  treating  it  as  6x6  elastic  dynamic 
stiffness  matrix  whose  elements  are  zero  except  for  the  (1,1),  (1,4),  (4,1)  and  (4,4)th 
elements  which  are  obtained  from  equation  (2.3.1).  For  the  special  case  in  which  fluid 
layers  form  consecutive  elements  it  will  be  necessary  to  set  to  zero  rows  and  columns 
corresponding  to  axial  and  tangential  displacements  at  the  common  interface,  the  diagonal 
elements  being  set  to  unity  to  allow  numerical  inversion  without  the  chore  of  condensing  out 
displacements  associated  with  these  null  rows  and  columns. 


2.5  Point  Force  Excitation 

A  point  force  applied  to  the  /th  interface  at  coordinate  (rQ,<|)j^  is  represented  by  the  3x1 
excitation  vector 

F(4»,z)  =  Fo5(z-zo)6(4>-«l>o)/ro,  (2-5.1) 

where  F(,=(F/,F^,f/)^  are  the  point  force  amplitudes.  The  spectral  representation  of  these 
forces  is 
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(2.5.2) 


=  (Fj/2nr</,)cos(mn(Zo+i<)/2I]exp(-i/i<t>Q), 


F,(rf^,n^)  =  (F//2itr^)sin[mii(Z(,+^)/2L]exp(-jrt<|>^ 

The  elements  of  the  force  vector  are  added  to  the  37  -2,  3y-l  and  3/th  elements  of  the 
excitation. 


The  spherically  spreading  far  field  sound  radiation,  in  spherical  polar  coordinates  (/?,0,4») , 
in  a  fluid  of  density  p  and  sound  speed  c,  is  obtained  from  cylindrical  baffle  theory  [1-3] 
as 


where 


0  A)  =  -tp<J^exp(t/tF)  ^  (-0l"lL4,(a)«,(^,n,iw)exp(t>Ki>) 
p^/c,  ,<|>;  ^kRs.m!d  ^  H'Jkbsmd) 


Am(cc)  =  exp(iaL)(mn/2L^)ll-(-l)"exp(-2iaL)]/(m^rt^/4L^-a^) 


(2.5.3) 

(2.5.4) 


is  the  inverse  Fourier  integral  transform,  divided  by  I ,  of  the  radial  displacement  function 
sm[mit(z*L)/2L]  which  is  assumed  to  be  zero  outside  the  interval  z-{-LJL) .  This  function 
must  be  evaluated  at  the  stationary  phase  wavenumber  o=Jfecos(0).  The  displacement 
is  the  third  element  in  the  solution  vector  U(«,m)  of  equation  (2.4.1),  ie  it  is  the 
spectral  amplitude  of  the  radial  displacement  of  the  outer  surface  of  the  layered  cylinder. 


2.6  Plane  Wave  Excitation 

A  plane  wave  arriving  from  angles  0=0^  and  <1>=<|),, 

p.  =  Pj,exp(-ifccsin0jCOS<|>j-tA:ysin0,sin<j>i-ifacos0j),  (2.6.1) 

is  expressed  in  cylindrical  polar  coordinates  as 

Pi  =  P,itvp{-ikzcos,^)Y,  (-O’"'exp(in(<|>-(|>.))J,„,(*rsin0.).  (2.6.2) 

A 

The  pressure  scattered  from  an  infinite  rigid  cylinder,  of  radius  b,  is 


Pr  =  -Potxpi-ikzcosQ)  5^  (-i)l*'exp(in(i|>-<|)j)) 

m 


y'l,|(««sin0.)/f,,|(*rsin0.) 


(2.6.3) 
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The  excitation  on  the  cylinder  surface  is  taken  as  the  blocked  pressure,  -(p,+p^ ,  evaluated 
at  r=b. 


Em;z)  = 


Pffixpi-ikzcosQ)  52 

II 


-2j(-0  l*'exp(  w(<|> 
«ibteine^',,,(Jk&sin0^) 


(2.6.4) 


The  Fourier  series  coefficients  of  this  excitation  confined  to  the  axial  coordinates  z=(-LX) , 
are  obtained  from  equation  (2.1.2),  using  the  axial  part  of  the  transform  alone,  as 


EibjiM  = 


=*cos0j)i(-Ol"lexp(-m4>p 

nikbsin6,lf'|,|(iy«in0.) 


(2.6.5) 


This  excitation  is  added  to  the  3rd  element  of  the  excitation  vector. 

The  pressure  scattered  to  the  far  field  in  the  exterior  fluid  is  the  sum  of  two  terms,  viz 

p(R,m  =  p„(i?,0,4>)  +  p,(/?.0.4>)  (2-6.6) 

where  is  the  pressure  scattered  by  a  rigid  cylinder  of  length  2Lj,  as  given  by  James  [3], 


p^(jR,0,<|>)  =  2iL4PoSin0;[sin(fcLj(cos0+cos0,)]/^j,(cos0+cos0j)]  x 

exp(ikR)  y  m(4>-4><)) 

nJl!sm0  B  /f'j^j(jfcbsin0) 


and  is  in  the  far  field  pressure  radiated  by  surface  vibrations  confined  to  z=(.-LJL).  This 
pressure  is  given  by  Equation  (2.5.3)  in  which  uj(,b/i^)  is  the  third  element  in  the  solution 
vector  U(n,m)  of  equation  (2.4.1).  Normally  the  length,  2L^,  for  calculating  the  "rigid" 
body  scattering  is  taken  as  the  length,  21,,  of  the  vibrating  cylinder.  The  pressure  p^  does 
not  obey  the  reciprocity  principle  as  0  and  0^  are  not  interchangeable.  This  is  a  defect  of 
cylindrical  baffle  theory  which  requires  that  caution  should  be  exercised  when  interpreting 
scattering  levels  as  0=180°  is  approached. 
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3.  Numerical  Examples 


3.1  General 

A  Fortran  program  (the  listing  is  not  given  here)  has  been  written  to  calculate  transfer 
functions  of  far  field  sound  level  at  selected  observation  angles  (6,<t>)  when  the  excitation  is 
a  point  force;  and  transfer  functions  of  both  monostatic  and  bistatic  target  strength  when  the 
excitation  is  a  plane  wave  at  arbitrary  incidence.  For  the  case  of  a  thin  steel  shell,  with  or 
without  a  fluid  coating,  excellent  numerical  agreement  has  been  obtained  with  results 
calculated  from  thin  shell  theories  (2-4],  but  these  comparisons  are  not  included  here. 
Transient  responses  are  readily  calculated  from  the  time-harmonic  transfer  functions  by  using 
an  auxiliary  program  for  which  the  time  variation  of  the  source  can  be  specified  by  a  library 
of  excitations;  see  the  report  by  James  [7]  for  details. 


3.2.  Radiation 

In  order  to  illustrate  an  application  of  the  program,  consider  a  plastic  cylinder  which  is 
reinforced  circumferentially  by  carbon  fibre,  and  covered  by  a  compliant  coating.  The 
constants  of  the  CFRP  cylinder  are  £^^=202.0,  v^j=0.32,  v^^=0.427, 

v,,^  =0.0117,  G^^=2.74,  ,  p=1520,  outer  radius  b=1.0,  inner  radius  a=0.95 

and  loss  factor  0.02.  The  compliant  coating  constants  are  £=0.0028,  v=0.40,  G=0.001 , 
p =1100.0,  outer  radius  6=1.05,  inner  radius  a =1.00  and  loss  factor  ti=0.1  .  All  values  are 
in  SI  units  with  the  Young’s  and  shear  moduli  being  in  GN/m^.  The  composite  cylinder  is 
immersed  in  water  for  which  p=1000  and  c=1500,  and  contains  a  light  fluid  whose  reaction 
on  the  shell  is  negligible.  The  excitation  is  a  unit  radial  point  force  located  at  the  centre  of 
the  cylinder  on  its  inner  surface,  at  (|>=0®  ;  the  far  field  sound  radiation  is  measured  at  beam 
aspect  <|)=0°,  0=90“. 


In  Figure  2(a)  is  shown  the  transfer  function  of  the  far  field  sound  level  in  the  absence  of  the 
cor  iliant  coating.  The  strong  resonances  can  be  identified  by  the  integer  wavenumbers 
where  m  is  the  number  of  half-axial  wavelengths  and  n  the  number  of  full 
circumferential  wavelengths,  though  this  is  not  done  here.  In  Figure  2(b)  is  shown  the 
transient  response  when  the  excitation  is  a  single  cycle  of  a  200  Hz  sine  wave.  The  small 
response  up  to  relative  time  r=10  ms  is  non-causal  due  to  the  abrupt  termination  of  the 
radiated  sound  transfer  function  at  5(X)  Hz  where  it  is  large;  formally  this  is  the  constraint 
that  a  function  cannot  be  both  time  limited  and  bandwidth  limited.  There  follows  a  rapid 
resonant  build-up  and  a  slowly  decaying  ringing  dominated  by  the  185  Hz  discrete  in  the 
transfer  function  spectrum. 

In  Figure  3(a)  is  shown  the  transfer  function  of  the  far  field  sound  level  when  the  CFRP 
cylinder  has  an  external  compliant  coating.  The  isolating  properties  and  large  damping  of 
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the  coating  have  reduced  overall  sound  levels  and  blunted  sharp  resonances.  In  Figure  3(b) 
is  shown  the  transient  response  when  the  excitation  is  a  single  cycle  of  a  200  Hz  sine  wave. 
The  reduction  in  amplitude  and  much  faster  decay  of  the  transient  far  field  pressure,  due  to 
the  coating,  is  most  evident  when  this  plot  is  compared  with  Figure  2(b). 
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4.  Conclusions 


Formulae  have  been  given  for  calculating  sound  radiation  from  and  sound  scattering  by  a 
laminated  composite  cylinder  whose  ends  are  simply  supported.  Numerical  results  of  both 
steady-state  and  transient  acoustics  have  demonstrated  an  application  of  the  model,  which  is 
believed  to  be  a  useful  tool  for  investigations  of  the  physical  mechanisms  that  control  the 
acoustics  of  finite  fibre  reinforced  shells. 
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Figure  1.  Layered  cylinder,  coordinate  systems  and  excitations.  The  cylinder  is  located 
between  rigid  baffles  which  are  not  shown. 
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Figure  2.  Beam  aspect  sound  radiation  from  a  circumfrentiaily  reinforced  CFRP  cylinder 
in  water:  (a)  transfer  function  of  a  unit  time-harmonic  point  force;  and  (b)  transient  sound 
radiation  when  the  excitation  is  a  single  cycle  of  a  200  Hz  force.  The  cylinder  is  uncoated. 
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Figure  3.  Beam  aspect  sound  radiation  from  a  circumferentially  reinforced  CFRP  cylinder 
in  water:  (a)  transfer  function  of  a  unit  time-harmonic  point  force;  and  (b)  transient  sound 
radiation  when  the  excitation  is  a  single  cycle  of  a  200  Hz  force.  The  cylinder  is  covered 
externally  by  a  soft,  highly  damped,  coating. 
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